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We summarize the status of calculations of the electroweak radiative 
corrections to W and Z boson production via the Drell-Yan mechanism 
at hadron colliders. To fully exploit the precision physics potential of the 
high-luminosity environment of the Fermilab Tevatron pp (Run II) and the 
CERN LHC pp colliders, it is crucial that the theoretical predictions are 
well under control. The envisioned precision physics program includes a 
precise measurement of the W boson mass and the (leptonic) weak mixing 
angle, as well as probing the Standard Model (SM) of electroweak inter- 
actions at the highest accessible center-of-mass energies. Some numerical 
results are presented. 
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1 Introduction 



The Standard Model of electroweak interactions (SM) so far withstood all ex- 
perimental challenges and is tested as a quantum field theory at the 0.1% level [[]]]. 
However, the mechanism of mass generation in the SM predicts the existence of a 
Higgs boson which, so far, has eluded direct observation. Direct searches at LEP2 
give a (preliminary) 95% confidence-level lower bound on the mass of the SM Higgs 
boson of Mh > 113.5 GeV 0. Indirect information on the mass of the Higgs boson 
can be extracted from the Mh dependence of radiative corrections to the W boson 
mass. With the present knowledge of the W boson and top quark masses, and the 
electromagnetic coupling constant, a(M§), the SM Higgs boson mass can be indi- 
rectly constrained to M H = 77tH GeV [|TJ by a global fit to all electroweak precision 
data. Future more precise measurements of the W boson and top quark masses are 
expected to considerably improve the present indirect bound on Mh'- with a precision 
of 30 MeV for the W boson mass, M\y, and 2 GeV for the top quark mass which are 
target values for Run II of the Tevatron [[J, M H can be predicted with an uncertainty 
of about 30%. In addition, the confrontation of a precisely measured W boson mass 
with the indirect SM prediction from a global fit to all electroweak precision data, 
M w = 80.385 ±0.022 GeV 0, will provide a stringent test of the SM. A detailed dis- 
cussion of the prospects for the precision measurement of Mw, and of the (leptonic) 
effective weak mixing angle, sin 2 at Run II and the LHC is given in Refs. [|3] 
and [|J, respectively. 

In order to measure Mw with high precision in a hadron collider environment it 
is necessary to fully control higher order QCD and electroweak radiative corrections 
to the W and Z production processes. The status of the QCD corrections to W and 
Z boson production at hadron colliders is reviewed in Refs. j5]||. Here we discuss 
the electroweak O(a) corrections to pp — > W ± — > V^Ui and pp — > 7*, Z — > l + l~ 
(I = e, /i) as presented in detail in Refs. MM and PJlQfl. 



2 Electroweak 0(ol) Corrections to pp — > W — ► I v 

The full electroweak 0(a) corrections to resonant W boson production in a general 
four-fermion process were calculated in Ref. [0] with special emphasis on obtaining a 
gauge invariant decomposition into a photonic and non-photonic part. It was shown 
that the cross section for resonant W boson production via the Drell-Yan mechanism 
at parton level, qiq v f f (7), can be written in the following form ||: 

^(o+D = ^(o) {1 + 2Ue{F^\s = M&) + F^{s = M^))\ 

+ £ [da® Ffafrt) + dot^ , (1) 

a— initial , / inal , 
inter f. 
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where the Born cross section, dM°\ is of Breit-Wigner form, and s and t are the 
usual Mandelstam variables in the parton center of mass frame. The (modified) weak 
corrections and the virtual and soft photon emission from the initial and final state 
fermions (as well as their interference) are described by the form factors F® eak and 
Fq ED , respectively. The IR finite contribution da%_+ 3 describes real photon radiation 
away from soft singularities. Mass singularities of the form ln(s/m^) arise when the 
photon is emitted collinear with a charged fermion and the resulting singularity is 
regularized by retaining a finite fermion mass (m/). F^^ al and da^^ al still include 
quark-mass singularities which need to be extracted and absorbed into the parton 
distribution functions (PDFs). The absorption of the quark-mass singularities into the 
PDFs can be done in complete analogy to gluon emission in QCD, thereby introducing 
a QED factorization scheme dependence. Explicit expressions for the W production 
cross section in the QED DIS and MS scheme are provided in Ref. ||. So far, in the 
extraction of the PDFs from data as well as in the PDF evolution, QED corrections are 
not taken into account. The latter result in a modified scale dependence of the PDFs, 
which is expected to have a negligible effect on the observable cross sections 0] . The 
numerical evaluation of the cross section is done with the parton level Monte Carlo 
program WGRAD HQ and results have been obtained for a variety of interesting W 
boson observables at the Tevatron || and the LHC 0. 

In the past, fits to the distribution of the transverse mass of the final-state lepton 
neutrino system, Mt(Iv), have provided the most accurate measurements of Mw [0]- 
Photonic initial state and initial-final state interference corrections were found to have 
only a small effect on the distribution, and weak corrections uniformly reduce the 
cross section by about 1%. However, final-state photon radiation significantly distorts 
the shape of the distribution, and thus considerably affects the extracted value of 
M\y- In the electron case, when taking into account realistic lepton identification re- 
quirements to simulate the detector acceptance, the electroweak radiative corrections 
are strongly reduced because electron and photon momenta are combined for small 
opening angles between the two particles. This eliminates the mass singular terms 
associated with final state radiation. The ratio of the full O (a 3 ) and lowest order dif- 
ferential cross section as a function of M^ilv) with and without lepton identification 
requirements taken into account is shown in Fig. [I]. 



A previous approximate calculation JL2] took only the real photonic corrections 



properly into account while the effect of soft and virtual virtual photonic corrections 
were estimated from the inclusive W — > lui^y) width. Weak corrections were ignored 
in Ref. | T2]| . Comparing the W mass shifts obtained using the calculations of Refs. [12 



and ||, one finds that the proper treatment of virtual and soft corrections and the 
inclusion of weak corrections induces an additional shift of O(10 MeV) in the extracted 
W boson mass. 



WGRAD is available from the authors. 
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Figure 1: The relative corrections to the Mt{Iv) distributions at the Tevatron when taking 
into account the full electroweak 0(a) corrections (from Ref. M). 



3 Electroweak 0(ql) Corrections to pp -> 7 *,Z -> l+l- 

Neutral-current Drell-Yan production is interesting for several reasons: 

1. Future precise measurements of the W boson mass at hadron colliders depend 
on a precise knowledge of the Z boson production process. When compared to 
the values measured at LEP, the measured Z boson mass and width help to 
determine the energy scale and resolution of the electromagnetic calorimeter. 

2. Ratios of W and Z boson observables may yield a more precise measurement of 
Mw than the traditional technique of fitting the Mt distribution @,[nj. 

3. The forward-backward asymmetry in the vicinity of the Z resonance can be 
used to measure the (leptonic) effective weak mixing angle |4]||. Studying the 
forward-backward asymmetry above the Z resonance probes the 7, Z interfer- 
ence at the highest available energies. 

4. Finally, at large di-lepton invariant masses, m(l + l~), deviations from the SM 
prediction could indicate the presence of new physics, such as new heavy gauge 
bosons Z' or extra spatial dimensions. 

It is therefore important to determine the electroweak corrections for this process. 
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The electroweak 0(a) corrections to neutral-current Drell-Yan processes naturally 
decompose into QED and weak contributions, i.e. they build gauge invariant subsets 
and thus can be discussed separately. The observable next-to-leading order (NLO) 
cross section is obtained by convoluting the parton cross section with the quark 
distribution functions q(x, Q 2 ) (s = XiX%S) fllOl 



da(S) =^ 1 rfx 1 c/x 2 g(x 1 ,g 2 )g(x 2 ,g 2 )[da(°+ 1 )(^t) + da Q ED(^ ED ,s,t)] , (2) 

where da ( - 0+1 - ) comprises the NLO cross section including weak corrections, and d<3"QED 
describes the QED part, i.e. virtual corrections and real photon emission off the 
quarks and charged leptons. The PDFs depend on the QCD renormalization and 
factorization scales which we choose to be equal; the common scale is denoted by 
Q 2 . The radiation of collinear photons off quarks requires the factorization of the 
arising mass singularities into the PDFs which introduces a dependence on the QED 
factorization scale, Hqed- The treatment of mass singularities is universal and thus 
the same as in the W case. The QED 0(a) corrections to pp — > j*,Z — > 
(I = e, fi) have been calculated and implemented in the parton level Monte Carlo 
program ZGRAD 0[] and their impact on the di-lepton invariant mass spectrum, the 
lepton transverse momentum distribution, and on the forward-backward asymmetry, 
Aps, has been studied. In addition, the prospects for a precision measurement of 
sin 2 #gjy extracted from Ap-Q at the Z resonance at the LHC were investigated. 

In Fig. H] we show the effect of the QED corrections on the invariant mass distri- 
bution of the final state lepton pair. Similar to the transverse mass distribution in 
the charged-current Drell-Yan process, final-state photon radiation strongly affects 
the shape of the m(l + l~) distribution. When lepton identification requirements are 
taken into account, the large contributions from mass singular logarithms largely 
cancel in the electron case. As in the charged-current Drell-Yan process, initial-final 
state interference is negligible, and the impact of initial-state radiation is small after 
factorizing the collinear singularities into the PDFs. The difference in the extracted 



Z boson mass when comparing the approximate calculation of Ref. |12j with the full 
calculation of the 0(a) QED corrections is of O(10 MeV). Since the detector response 
is calibrated using Z boson observables, the shift in the Z boson mass is expected to 
slightly modify the W mass extracted from experiment. 

For precision physics away from the Z resonance, the (non- universal) weak cor- 
rections must also be included. These corrections become important at large values 
of the di-lepton invariant mass due to the presence of large Sudakov-like electroweak 
logarithms of the form \n(m(l + l~)/My), V = W, Z, which eventually may be re- 
summed [[TH]. A calculation of the non-universal weak corrections in pp — ► 7*, Z — ► 



/ + / is currently in progress |T(J . In the implementation of the weak corrections we 



2 ZGRAD is available from the authors. 
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Figure 2: The relative corrections to the m(e + e ) and m(^ + ^ ) distributions in Drell-Yan 
production at the Tevatron due to the 0(a) QED corrections (from Ref. ||). 



closely follow Ref. |Tj|, in particular for the treatment of higher-order corrections, 



which are important for a precise description of the Z resonance. 

The NLO parton differential cross section, including weak O(o) and leading 0(a 2 ) 
corrections, which enters eq. (H) is of the form 



da^ = dP 2f 1 £ |4 0+1 ) + 4° +1) | 2 (M) + da hox (s,t) , (3) 

where the sum is taken over the spin and color degrees of freedom, and dP2f denotes 
the two-particle phase space. d<5"box describes the contribution of the box diagrams 
involving two massive gauge bosons. The matrix elements A^ 1 ^ comprise the Born 
matrix elements, the 7, Z, ^Z self energy insertions including a leading-log resumma- 
tion of the terms involving the light fermions, and the one-loop vertex corrections. 
While A^ 1 ^ can be expressed in terms of effective vector and axial-vector couplings, 
the box contribution dow cannot be absorbed in effective couplings. However, in 
the Z resonance region the box diagrams can be neglected and the NLO cross section 
d<j( 0+1 ) of eq. (§) has a Born-like structure. The leading universal electroweak correc- 
tions, i.e. the running of the electromagnetic charge and corrections connected to Ap, 
can be included in form of an effective Born approximation (EBA). Comparing results 
of the calculation which includes the full 0(a) corrections with those obtained using 
the EBA together with the pure QED corrections reveals the effects of the genuine 
non-universal electroweak corrections such as box diagrams. 

The weak corrections to neutral-current Drell-Yan processes as described above 
are currently being implemented in the parton level MC program ZGRAD2 fL0|. A 



detailed numerical discussion of the effects of the electroweak 0(a) corrections on 



distributions in pp — > 7*, Z — > / + Z (7), I — e, /J, at the Tevatron and the LHC will be 
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Figure 3: The relative corrections to the m(^ + /i~) distribution a) at the Tevatron and b) 
at the LHC when taking into account the universal corrections entering the EBA and QED 
corrections only (solid line), and when the full O(a) electroweak corrections are included in 
the calculation (dashed line). 



given in Ref . [|T(J] . Here we present some selected preliminary results for the di-lepton 
invariant mass distribution and the forward-backward asymmetry. 

In Fig. |3| we show the yU + /i~ invariant mass distribution including the full 0(a) 
corrections normalized to the differential cross section in the EBA for large di-lepton 
invariant masses at the Tevatron and the LHC. Separation cuts and lepton iden- 
tification requirements to simulate the detector acceptance as described in Ref. 
(Tevatron) and Ref. [|J] (LHC) are taken into account in Fig. [| For comparison the 
relative corrections including the QED corrections only are also shown. As expected 
from the presence of large electroweak Sudakov-like logarithms, the weak corrections 
strongly increase in magnitude with increasing m(/i + /i~), reaching about 10% at 
m(/z + yU~) = 1 TeV. Both, the QED and the genuine weak corrections reduce the 
differential cross section. Qualitatively similar results are obtained in the e + e~ case. 

In Fig. U we show how the purely weak corrections affect the forward backward 
asymmetry at the LHC[]. To illustrate the effect of the non-universal weak corrections, 
we plot the difference of the forward backward asymmetry including the full 0(a) 
corrections, and the asymmetry which only takes into account QED corrections and 
the universal corrections which are included in the EBA. A genuine non-universal 
electroweak effect can be observed in the vicinity of m(l + l~) = M\y and 2 M\y, which 
is due to threshold effects in the box diagrams involving two W bosons. Results 
qualitatively similar to those shown in Fig. [| are also obtained for the Tevatron. 

3 For a definition of ^4fb at the LHC, see Ref. [||. 
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Figure 4: The forward-backward asymmetry including NLO electroweak corrections at the 
LHC, imposing the cuts and lepton identification requirements of Ref. |T6| ]. The EBA and 
QED contribution have been subtracted (preliminary results). 



The forward backward asymmetry at the LHC is very sensitive to the rapidity 
coverage of the leptons assumed. In Fig. f|, we have used the lepton rapidity coverages 
foreseen for the ATLAS detector j^Io"]. For muon pairs, both muons are required to 



have rapidity \y(/^)\ < 2.5. For e + e~ pairs, the leptons are required to have \yi(e) \ < 
4.9, with one of them having to fulfill the more stringent requirement ||/f(e)| < 2.5. 
In addition, the lepton pair rapidity has to be \y(ll)\ > 1 for both electrons and 
muons in the final state. This cut substantially increases the magnitude of Ap-Q at 



the LHC UTTf 



It is interesting to check whether the threshold effect at m(l + l ) = 2 Mw will 
be observable. In the electron case, the expected statistical uncertainty in A^b for 
m(e + e~) = 2 M w ± 5 GeV and 100 fb" 1 at the LHC is about (3 - 4) x 10" 3 per 
experiment. The size of the non-universal electroweak corrections in the region are of 
the order of 10~ 3 . In a realistic calculation, contributions from W + W~ — > l + vil~Vi, 
ZZ — > l + l~Vv and tt production to the forward backward asymmetry need to be taken 
into account, which could well be of the same order as the genuine weak corrections. 
It will thus be difficult to observe a clear signal of the threshold effects originating 
from the box diagrams involving two W bosons in A-ps at the LHC. On the other 
hand, given the expected statistical precision, the genuine weak corrections cannot 
be neglected when comparing data with the SM prediction. 
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4 Conclusions 



Our results show that, for the precision obtained in previous Tevatron runs, the 
existing calculations for W and Z boson production are sufficient. However, for fu- 
ture precision measurements the full electroweak 0{a) corrections and probably also 
multiple photon radiation effects should be taken into account. The inclusion of the 
non-resonant contributions to W production in WGRAD is in progress fl4| (see also 
Ref. |19|]). As a first step towards a calculation of the 0(a 2 ) QED corrections, the 



effects of two-photon radiation in W and Z boson production at hadron colliders have 
been computed in Ref. pO 
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